Stable two-dimensional dispersion-managed soliton 
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The existence of a dispersion-managed soliton in two-dimensional nonlinear Schrodinger equation 
with periodically varying dispersion has been explored. The averaged equations for the soliton width 
and chirp are obtained which successfully describe the long time evolution of the soliton. The slow 
dynamics of the soliton around the fixed points for the width and chirp are investigated and the 
corresponding frequencies are calculated. Analytical predictions are confirmed by direct PDE and 
ODE simulations. Application to a Bose-Einstein condensate in optical lattice is discussed. The 
existence of a dispersion-managed matter- wave soliton in such system is shown. 
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I. INTRODUCTION 

Nonlinear wave propagation in media with periodi- 
cally varying dispersion is attracting a huge interest over 
the recent years. A prominent example is a dispersion- 
managed (DM) optical soliton, which is considered to be- 
come the major concept in future soliton-based commu- 
nication systems. It was shown theoretically and experi- 
mentally that the strong DM regime provides the undis- 
turbed propagation of pulses over very long distances. 
DM solitons are robust to the Gordon-Haus timing jit- 
ter, which makes them favorable against the standard 
solitons [1,2]. 

Mathematically this type of problems are described by 
the one dimensional (ID) nonlinear Schrodinger (NLS) 
equation with periodic dispersion - a nonlinear analogue 
of the Mathieu equation. The corresponding linear equa- 
tion exhibits a rich variety of stability and instability 
zones for the parameters. The existence of a DM soli- 
ton is one of the nontrivial consequences of the stable 
diagram for the periodic NLS equation. 

Although well studied in the ID case, the two and three 
dimensional extensions of this problem are far less ex- 
plored. The major difference here is that, contrary to 
the ID case, the NLS equation in two and three dimen- 
sions is unstable against collapse. In particular, for the 
two dimensional (2D) case the collapse occurs if the ini- 
tial power exceeds some critical value, i.e. if E > E cr . 
Recently it has been demonstrated that the nonlinear- 
ity management can prevent the collapse of solitons in 
2D Kerr type optical media [3,4], as well as in 2D Bose- 
Einstein condensates [5,6]. From these one can reason- 
ably expect that the dispersion-management can play 
balancing role also in the 2D case, and the stable 2D 
DM soliton can exist. Such a possibility has recently 
been considered in Ref. [7] by construction of the ground 
state for the periodic 2D NLS equation based on the av- 



eraged variational principle and the techniques of inte- 
gral inequalities - i.e. the proof of the existence theorem 
for DM soliton was presented. Analytical and numer- 
ical treatment of the problem, however, has not been 
addressed so far. 

The purpose of this Communication is to derive ana- 
lytical expressions for the parameters of a 2D dispersion- 
managed soliton and to study the conditions for their 
stability. To this regard, we use a time-dependent vari- 
ational approach (VA) to derive a set of ODEs for the 
soliton parameters. The stability of the DM soliton is 
then inferred from the stability of fixed points of the VA 
equations. 

The field dynamics is governed by the following 2D 
NLS equation 

iu t + d(t)Au + \u\ 2 u = 0, (1) 

where d(t) = do + d±(t) represents a time periodically 
varying dispersion coefficient. In the strong DM regime 
it is assumed that d(t) ~ (1 / e)d(t / e) , e <C 1 and the dis- 
persion averaged over the period is < d(t) >= do (in 
this case do > corresponds to a negative dispersion and 
d < to a positive one) . 

The equation (1) can be associated with two main 
physical problems: (i) beam propagation in 2D waveg- 
uide arrays with periodically variable coupling between 
waveguides [8,9]; (ii) nonlinear matter-waves of Bose- 
Einstein condensates in 2D optical lattices. 

In case (i) the model equations for a 2D nonlinear fiber 
array are given by [10] 

ilpn,z + K(z)A 2 ^ n + Uj"lpn,tt + XlVU Vn = 0, (2) 

where ip n is the envelope of electric field in the n-th fiber, 
A2 is the finite second difference for 2D, k(z) is the vari- 
able along z coupling coefficient [8,9], u>" is the group- 
velocity dispersion, x is the coefficient of nonlinearity. 
For long wavelength pulses the group-velocity dispersion 
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lu" can be neglected. Introducing the dimensionless vari- 
ables kz = t,ip n = \j2n/xu n and considering the field 
distribution to be broad in the transverse direction (> 7 
sites), one arrives to Eq. (1) with time and space in- 
terchanged and with d(t) describing a varying diffraction 
along the longitudinal direction. Notice that although 
the intrinsic discreteness of the array may arrest the col- 
lapse of a 2D NLS wave, it does not necessarily stabilize 
the pulse against decay. In the following we show that 
this can be done employing dispersion (diffraction) man- 
agement by means of which a stable 2D soliton can be 
created before the strong shrinking of the wave occurs. 

A similar situation arises in case (ii) for a Bose-Einstein 
condensates (BEC) confined in a 2D optical lattice. In 
this case dynamics of the condensate is governed by the 
Gross-Pitaevskii (GP) equation 

h 2 

im t = -— A* + V{x, y; t)V + 52 d|*| 2 *, (3) 
2m 

where g 2 D = 93d/ V2^a z ), g 3D = 4irh 2 a s /m,a z = 
(H/mu} z ) 1/2 and with V(x,y) = V (t)(cos 2 (k x) + 
cos 2 (fcoy)) denoting an optical lattice with the ampli- 
tude periodically varying in time. Spatiotemporal wave 
collapse in the framework of a similar equation (when 
the potential is periodic in one direction V(x, y) — 
Vocos(fcoa;)) was considered in [11], where analytical ex- 
pression for the upward shift of collapse criterion was 
derived for potentials rapidly oscillating in space (large 
fco)- By adopting an effective mass description one can 
show that the 2D GP equation [12] can be reduced to 
the DM NLS equation (1). The effectiveness of DM ap- 
plied to quasi-lD atomic matter-waves was experimen- 
tally demonstrated in Ref. [13]. 

For analytical considerations it is convenient to refer 
to the axially symmetric case for which A = d 2 /dr 2 + 
(l/r)(d/dr), and apply the harmonic modulation for 
dispersion-management: d(t) = d + d\ sin(Qt). We re- 
mark that although in the present Communication we do 
not consider the case of two-step dispersion-management: 
d(t) =d+, if t + + nt p >t> nt p , and d(t) = d—, if (n + 
l)t p > t > ntp + t_|_, where t p = t+ + t_, n = 0, 1, 2... , 
this approach can also be effectively used for the creation 
of stable 2D DM solitons. 

Our analysis of the pulse dynamics under dispersion- 
management is based on the variational approach [2,14], 
according to which a space averaged Lagrangian L = 
J Ldf is constructed starting from a suitable ansatz for 
the soliton profile. In the following we shall calculate L 
by using the following Gaussian ansatz 

r 2 b(t)r 2 

u(r,t) = A(t)exp(-^+i-^— +i#t)), (4) 

where A, a, b, <j> denote the amplitude, width, chirp and 
linear phase of the soliton, respectively. The equations 
for the soliton parameters are then derived from the 
Euler-Lagrange equations for L as 



<H = 2d(t)0, &= 2d{t) a s E > (5) 
where (5 — ab, and E = J °° \u\ 2 rdr is the energy. 



II. SYSTEM OF AVERAGED VARIATIONAL 
EQUATIONS 

Let us consider the evolution of a pulse (a beam or a 
soliton matter-wave, depending on the physical system 
in consideration) using the division on the fast and slow 
time scales [15-17]. The width and chirp of the pulse 
are then represented as a{t) = a + oi, (3(t) = /3 + 
where a, (3 are slowly varying functions on the scale 1/e 
and ai,p\ are rapidly varying functions. The solutions 
for ai,p\ are 

4eMi . 2Qd!P 
a i = — 9 — ; — 7^o\ sm(s2t) - , i2 , ^ 2 cos ("*)> ( 6 ) 



/3i = 



a 3 (co 2 + Q 2 ) 
6adi/3 



d 4 (u; 2 + n 2 ) 



sm(Ctt) 



2difl 



o?{oj 2 + n 2 ) 



cos(m), (7) 



where uiq = —6a/a 4 , a = 2d — E. Note that a < for 
over-critical energy for collapse E > E cr = 2, at do = 1 
given by the VA. The exact value, corresponding to the 
so called "Townes soliton" is E cr — 1.862 [18]. Consid- 
ering the limit of high frequencies SI 2 3> uj 2 <~ 1 , for the 
averaged parameters of the system we finally get 



d t = 2f3{d 



I2d 2 d n 

n 2 a 7 
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n 2 a 4 
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n 2 d 5 



), 



(8) 
(9) 



This system has the Hamiltonian structure with the 
Hamiltonian given by 



a 2A 2 d -,, , 



from which the equations of motion follows as dt = 
dH/d/3, (3t — —dH/dd. From this Hamiltonian one can 
also see that the mechanism for collapse suppression orig- 
inates from the repulsive potential near the small values 
of width ~ 1/a 6 , which counteracts to the attractive force 
induced by the nonlincarity <~ 1/a 2 . The exact balance 
between these forces gives rise to a stable state. This 
state is oscillatory with the frequency which will be de- 
fined later. The stabilization mechanism of a 2D NLSE 
soliton is similar to that of the inverted pendulum with 
oscillating pivot point [19]. We should note that the av- 
eraged dynamics is not potential - a velocity dependent 
term appears in the interaction potential (see 4th term in 
(10)). Although this term doesn't contribute to the fixed 
point, it is important for the description of oscillatory 
dynamics of 2D DM solitons. 

The system (8), (9) has the fixed points 
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= 0, a c = (_l^!)V4. 



(11) 



employed, which also imitates the infinite domain condi- 
tion. The dispersion map was supposed to have parame- 
ters do = 1, d\ = 3.5, O = 50. 
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FIG. 1. Phase portrait of the variational equations (5) with 
parameters do = 1, di = 3.5, Q. = 50, E — N/2n = 2.3034. 

Note that A is proportional to the strength of the dis- 
persion map D = 2ndi/n, therefore a c <~ \J~D in analogy 
with the estimate for a DM soliton in ID case. There ex- 
ists one solution with a stationary width for the anoma- 
lous residual dispersion e? > 0, E > 2d . This is con- 
firmed by the phase portrait (Fig.l) of the variational 
system (5). 

Let us analyze the stability of fixed points for the 
anomalous residual dispersion do > 0. We assume 
a = a c + ea\,(3 = e(3\. Substituting into Eq.(8) and 
Eq.(9), and collecting terms of order e we find 



ai,t 



o a 2 

(2d + -^)/3i =M(3 1 , 



a 

,3<r 84A% 
Pi,t = -(-4 + — — )ai 



-Sai. 



(12) 
(13) 



The oscillations of the width and chirp near the fixed 
points are stable if MS > 0, which is always satisfied 
for d > 0, E > 2d . The frequency of secondary slow 
oscillations of a 2D DM soliton is proportional to >/MS. 



III. NUMERICAL SIMULATIONS 

To avoid the singularity at r — we consider the 
problem in Cartesian coordinates A = d 2 + d 2 , and 
r 2 = x 2 + y 2 . Then numerical simulations can be per- 
formed by two-dimensional fast Fourier transform [20]. 
The results are produced using a 2D grid of 256 x 256 
points over the domain x,y e [—6.4,6.4] and the time 
step 5t = 0.001. To prevent the back-action of a small 
amount of linear waves, resulted from the periodic per- 
turbation, the absorption on the domain boundaries is 



FIG. 2. Evolution of a 2D DM soliton according to numer- 
ical solution of the equation (1). The wave function is nor- 
malized to N = 2nEo with Eq = 2.3034, and the dispersion 
map is do = 1, di = 3.5, = 50. 



This choice of parameters corresponds to moderate 
dispersion- management (D ~ 0.45). The axial section 
profile of the wave function \u(r, t)\ 2 as obtained by di- 
rect numerical solution of the PDE (1) is presented in 
Fig. 2. As can be seen, rather stable quasi-periodic dy- 
namics is realized for a selected parameter settings. Note 
that would the periodic modulation of the dispersion had 
not been applied, the initial waveform would have col- 
lapsed within t <~ 3. The dispersion-management sta- 
bilizes the pulse against the collapse or decay, provid- 
ing undisturbed propagation over very long distances. 
The agreement between the predictions of the variational 
equations (5) for the width of a 2D DM soliton and the 
corresponding result from the full PDE simulations is re- 
ported in Fig. 3. As can be observed from this figure, 
the width of a 2D DM soliton performs quasi-periodic 
motion with the average width of a ~ 0.8 according to 
variational equations, while the PDE simulation yields 
a ~ 0.7. The fixed point for the above set of parameter 
values, according to eq.(ll) is a c = 0.6635 (see Fig.l). 
The frequencies of slow dynamics given by the VA equa- 
tions and PDE arc also in well agreement (Fig. 3). The es- 
timate for the frequency of slow oscillations from Eq.(12) 
yields uj a = >/MS = 3.5, therefore the period is T a = 1.9. 
The direct gauge from the Fig. 3 shows that T a ~ 2.2, in 
reasonable agreement with the above VA estimate. 

For Bosc-Einstcin condensates in a 2D optical lattice 
the dispersion coefficient can be expressed as d(t) = 
m/m*(t) in the effective mass formalism [12]. The ef- 
fective mass m* substantially differs from the true mass 
to (becoming even negative) and can be varied by chang- 
ing the parameters of the periodic potential, or inducing 
the transitions between energy bands. 
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FIG. 3. Stable quasi-periodic dynamics of the width of 
a 2D DM soliton. Solid line - variational equations (5) 
solved for E = N/2-k — 2.3034, and the initial conditions 
a(0) = 1, (3(0) = 0. Dashed line - full PDE simulations of the 
equation (1). 

For example, transitions between the 1st and 2nd 
bands (at the band edges) in the optical lattice of 
strength Vq = 2AE rec ( where E rec = fi'k^jlra is the 
recoil energy, ko = 2ir/\o, Xq is the laser wavelength), 
leads to variation of the dispersion coefficient in the range 
d(t) = (—2.5 4- 4.5) as considered above. 

IV. CONCLUSION 

In conclusion, we have demonstrated the possibility to 
stabilize the 2D soliton with over-critical energy (E > 
E cr ) by applying the dispersion-management. The de- 
veloped theory based on the variational approximation 
successfully describes the long term evolution of a 2D 
DM soliton, which is confirmed by direct PDE simula- 
tions. We discussed the possible experimental realiza- 
tion of a stable 2D dispersion-managed soliton in Bosc- 
Einstcn condensates confined to optical lattices. 
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